A Search for WIMPs with the First Five-Tower Data from CDMS 
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We report first results from the Cryogenic Dark Matter Search (CDMS II) experiment running 
with its full complement of 30 cryogenic particle detectors at the Soudan Underground Laboratory. 
This report is based on the analysis of data acquired between October 2006 and July 2007 from 15 
Ge detectors (3.75kg), giving an effective exposure of 121.3 kg-d (averaged over recoil energies 10- 
100 keV, weighted for a weakly interacting massive particle (WIMP) mass of eOGeV/c^). A blind 
analysis, incorporating improved techniques for event reconstruction and data quality monitoring, 
resulted in zero observed events. This analysis sets an upper limit on the WIMP-nucleon spin- 
independent cross section of 6.6x lO"*** cm^ (4.6x 10"'*'* cm'^ when combined with previous CDMS 
Soudan data) at the 90% confidence level for a WIMP mass of 60GeV/c^. By providing the best 
sensitivity for dark matter WIMPs with masses above 42 GeV/c^ , this work significantly restricts 
the parameter space for some of the favored supersymmetric models. 

PACS numbers: 14.80.Ly, 95.35.+d, 95.30.Cq, 95.30.-k, 85.25.0j, 29.40.Wk 
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Cosmological observations J] imply the existence of 
non-baryonic dark matter that drives structure forma- 
tion on large scales and dominates galactic and extra- 
galactic dynamics. Weakly Interacting Massive Particles 
(WIMPs) 3, with masses between a few tens of GeV/c^ 
and a few TeV/c^, form a generic class of dark matter 
candidates, motivated 0, 13| both by the measured value 
of the cosmological density and by the need to stabilize 
the standard model of particle physics at the weak scale. 

WIMPs should be distributed in the halo surrounding 
the Milky Way and scatter in terrestrial particle detec- 
tors [5j, 13]. Their coherent scattering on nuclei should 
lead to a roughly exponential energy-transfer spectrum 
with a mean recoil energy in the tens of keV 0,13| • The 
event rate is expected to be below 0.1 event per kilogram 
of target per day, much smaller than radioactivity rates in 



most materials. A number of technologies, most based on 
the identification of nuclear recoils among the electron re- 
coils produced by gammas and betas from radioactivity, 
are starting to reach this sensitivity level, corresponding 
to a spin-independent WIMP-nucleon scattering cross- 
section on the order of 10"^'^ cm^. Such "direct" searches 
for WIMP elastic scattering are complementary to "in- 
direct" searches for their annihilation products in our 
galaxy and to searches for new physics (e.g. supersym- 
metry or additional dimensions) at particle colliders [8] . 

The Cryogenic Dark Matter Search (CDMS II) oper- 
ates a total of 19 Ge (250 g each) and 11 Si (100 g each) 
solid-state detectors at ~40mK in the Soudan Under- 
ground Laboratory @,[l3|- Each detector is a disk, 7.6 cm 
in diameter and 1 cm thick. Ionizing radiation produces 
electrons and holes together with phonons. These charge 
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carriers are drifted by a small electric field (3 V/cm) and 
collected on two concentric electrodes on one flat face. 
Athermal phonons are collected using four circuits of su- 
perconducting thin-films, each circuit covering a quad- 
rant on the other flat face. The ratio of ionization to 
phonon recoil energy ("ionization yield") allows us to 
discriminate nuclear from electron recoils with a rejec- 
tion factor of > 10^. Electron recoils within ~10/xm of 
the detector surface suffer from a suppressed ionization 
signal. The resulting reduction in ionization yield can 
be sufficient to misclassify a surface electron recoil as a 
nuclear recoil. Signal timing provides effective discrim- 
ination against these events, improving our overall re- 
jection of electron recoils to > 10^. To reduce the ex- 
ternal gamma and neutron backgrounds, the experimen- 
tal setup [T^ also includes passive Pb and polyethylene 
shielding, which is surrounded by an active scintillator 
veto to detect cosmogenic muons or showers which could 
produce events in our detectors. 

We report on data from two periods (designated as 
runs 123 and 124) between October 2006 and July 2007. 
Improvements made since our previous publications 11 1 
include deployment of 18 additional detectors, increased 
exposure, greater cryogenic stability, faster data acqui- 
sition, enhanced monitoring and control of data quality, 
and improved analysis techniques. In this analysis we 
consider only the Ge detectors for WIMP search. Of 
the 19 Ge detectors, three suffering reduced performance 
from readout failures and one with relatively poor energy 
resolution have been left out of the present report. The 
remaining 15 Ge detectors (3.75 kg) were used for the 
run 123 analysis. Eight of these detectors were excluded 
from this WIMP search analysis during the shorter run 
124 due to variations in performance between the two 
runs. The data from these excluded Ge detectors and 
the Si detectors will be analyzed for later publication. 

Extensive calibrations with gamma (^'^^Ba) and neu- 
tron (^^^Cf) sources were used to determine data selec- 
tion criteria ( "cuts" ) that define the signal region, mon- 
itor detector stability, and characterize detector perfor- 
mance. Calibration data taken regularly with ^^'^Ba gen- 
erated over 28 million electron-recoil events between 10- 
100 keV, exceeding by a factor of thirty the number of 
comparable events in the WIMP-search data. Alternat- 
ing events from the ^^'^Ba data were separated into two 
statistically independent samples to allow unbiased char- 
acterization of cut performance. Over 600,000 events 
were recorded using the ^^^Cf source during five separate 
periods throughout the runs, including more than 10^ 
nuclear recoils used to characterize WIMP acceptance. 

Both calibration and WIMP search data were used to 
study detector stability and identify periods of poor per- 
formance. A standard set of one- and two-dimensional 
event parameter distributions were identified and data 
sets with significant deviations were discarded. A cut 
was also developed to remove periods of temporary poor 



ionization collection, during which bulk-electron recoils 
may leak into the signal region. Finally, we remove events 
from the outer part of each detector by a fiducial volume 
cut based on the partitioning of energy between the two 
concentric charge electrodes. 

The phonon energy scale and the timing of the phonon 
signals vary slightly depending on the position and en- 
ergy of each event. In order to maintain effective surface 
event rejection we compensate for these variations using 
an empirical look-up table based on our electron-recoil 
calibrations. The present analysis incorporates energy 
dependencies into this look-up table alongside position 
dependencies for the first time, enabling improved sur- 
face event discrimination. 

Event reconstruction at large radius remains imper- 
fect due to degeneracies in the phonon position quantities 
that inform the look-up table. A small number of high- 
radius events suffered from miscalibration due to these 
degeneracies. We developed a cut on the position-related 
phonon quantities from calibration data to remove events 
in problematic regions of the look-up table. 

We require that a candidate dark matter event deposit 
significant energy (> 4ct above mean noise) in one and 
only one detector ("single scatter event"), since WIMPs 
will not interact more than once in our apparatus. All 30 
detectors contributed to active vetoing of multiple scatter 
events at all times. We further require that a WIMP 
candidate show no significant activity in the surrounding 
scintillator veto shield during a 200 /us window around 
the trigger. Candidates must also lie within the 2a region 
of the nuclear-recoil distribution in ionization yield as 
determined by neutron calibration. 

For each event, we measure the risetime of the largest 
phonon pulse, and also its delay relative to the faster ion- 
ization signal. A cut based on the sum of the risetime and 
delay provides good rejection of surface electron-recoil 
events while retaining reasonable acceptance of nuclear 
recoils. Figure [1] shows an example of the distribution 
of this "timing parameter" in calibration data (gammas, 
gamma-induced surface events and neutron-induced nu- 
clear recoils). To effectively remove surface events we 
require that candidate dark matter events exceed a min- 
imum value for the timing parameter ("timing cut"), de- 
termined individually for each detector by setting an al- 
lowed maximum passage fraction for surface events in a 
subset of the ^'^^Ba calibration data. We also require that 
WIMP candidates be consistent with the nuclear-recoil 
event distribution (i.e. the difference between delay and 
risetime is less than a 4(t deviation from the neutron dis- 
tribution mean) . The performance of this cut is superior 
to that of earlier analyses due to improvements to the 
look-up table. 

The acceptance of our analysis cuts for single-scatter 
nuclear recoils was measured as a function of energy 
based on both neutron calibration and WIMP search 
data. Most cuts have very little effect on our acceptance 
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FIG. 1: Ionization yield versus timing parameter (see text) for 
calibration data in one of our Ge detectors. The yield is nor- 
malized to unity for typical bulk-electron recoils (dots; from 
^^^Ba gamma rays) . Low-yield ^'^'^Ba events (-I-) , attributed to 
surface electron recoils, have small timing parameter values, 
allowing discrimination from neutron-induced nuclear recoils 
from ^^^Cf (o), which show a wide range of timing parameter 
values. The vertical dashed line indicates the minimum tim- 
ing parameter allowed for candidate dark matter events in this 
detector, and the box shows the approximate signal region, 
which is in fact weakly energy dependent. (Color online.) 

of true nuclear recoils, with the ionization-based fidu- 
cial volume and phonon-timing cuts imposing the high- 
est costs in signal acceptance, both measured on neutron 
calibration data, as shown in Figured The exposure of 
this analysis is 397.8kg-days before and 121.3kg-days af- 
ter these cuts (averaged over recoil energies 10-100 keV, 
weighted for a WIMP mass of 60 GeV/c^). 

To avoid bias, we performed a blind analysis. An event 
mask was defined during initial data reduction to re- 
move events in and near the signal region from WIMP 
search data sets while developing the analysis. This mask 
was based on primary quantities not subject to refine- 
ment during the analysis process, keeping the event selec- 
tion constant throughout the analysis process described 
above. After WIMP selection criteria were finalized, the 
masking was relaxed to cover only the actual signal re- 
gion to aid in background estimation. 

Surface events mainly occur due to radioactive contam- 
ination on detector surfaces, or as a result of external 
gamma ray interactions releasing low-energy electrons 
from surfaces near the detectors. A correlation analy- 
sis between alpha-decay and surface-event rates provides 
evidence that ^^°Pb (a daughter-product of ^^^Rn) is a 
major component of our surface event background [l^ . 
Surface events generated in situ at Soudan, either from 
calibration with a ^^'^Ba source or naturally present in 
the WIMP search data, were studied to understand the 
surface event leakage into the signal region. We estimate 
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FIG. 2: Nuclear-recoil acceptance efficiency for event-specific 
cuts (i.e. excluding discarded data periods) as a function of 
recoil energy, averaged over all detectors used in the current 
analysis, weighted by their individual livetimes. The four 
curves represent the cumulative efficiencies at various stages 
during the analysis, culminating with the final efficiency (bot- 
tom) used to generate Figure |4l 

the surface event leakage based on the observed num- 
bers of single- and multiple-scatter events in each detec- 
tor within and surrounding the 2a nuclear-recoil region. 
The expected background due to surface interactions in 
this WIMP search analysis is 0.6 ± 0.5 events. 

Neutrons induced by radioactive processes or by 
cosmic-ray muons interacting near the apparatus can 
generate nuclear-recoil events that cannot be distin- 
guished from possible dark matter interactions on an 
event-by-event basis. Monte Carlo simulations of the 
cosmic-ray muons and subsequent neutron production 
and trans por t have been conducted with FLUKA [l3 |. 



MCNPX [111 and GEANT4 [Ij to estimate this cos- 
mogenic neutron background. Normalizing the results 
to the observed veto-coincident multiple-scatter nuclear- 
recoil rate leads to a conservative upper limit on this 
background of <0.1 events in our WIMP-search data. 

Additional Monte Carlo simulations of neutrons in- 
duced by nuclear decay were based on gamma-ray mea- 
surements of daughter products of U and Th in the ma- 
terials of our experimental setup and the assumption of 
secular equilibrium. The respective background estimate 
is <0.1 event, dominated by the deduced upper limit of 
U in the Pb shield. Direct measurements of U in Pb [l^ 
from the same source as the Pb used in our shield suggest 
a considerably lower contamination. 

After all analysis cuts were finalized and leakage esti- 
mation schemes selected, we unmasked the WIMP search 
signal region on February 4, 2008. No event was observed 
within the signal region. Figure [3] is a compilation of the 
low-yield events observed in all detectors used in this 
analysis. The upper panel shows the ionization yield dis- 
tribution versus energy for single-scatter events passing 
all data selection cuts except the timing cut. The four 
events passing the timing cut shown in the lower panel 
are outside the 2a nuclear-recoil region. 
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FIG. 3: Top: Ionization yield versus recoil energy in all de- 
tectors included in this analysis for events passing all cuts 
except the ionization yield and timing cuts. The signal re- 
gion between 10 and 100 keV recoil energies was defined using 
neutron calibration data and is indicated by the curved lines. 
Bulk-electron recoils have yield near unity and are above the 
vertical scale limits. Bottom: Same, but after applying the 
timing cut. No events are found within the signal region. 

Figure [?] shows the Poisson 90% C.L. upper hmit on 
the spin-independent WIMP-nucleon cross section de- 
rived from this data set (upper solid curve), based on 
standard assumptions about the galactic halo [7]. The 
minimum hes at 6.6xl0^'''*cm^ for a 60GeV/c^ WIMP. 

Our previous data from Soudan [lol [lH have been re- 
analyzed J_3] yielding a slight improvement in sensitiv- 
ity over our previous publications (upper curve in Fig- 
ure U]). A combined limit from all Soudan data (lower 
solid curve in Figure S]), using Yellin's Optimum Interval 
method [1S\ to account for observed events, gives an up- 
per limit of 4.6 x 10"'*'' cm^ at 90% C.L. for a WIMP mass 
of 60GeV/c^, a factor of ^3 stricter than our previously 
published limit. 

We also analyzed our data in terms of spin-dependent 
WIMP-nucleon interactions. Under the assumption of 
spin-dependent coupling to neutrons alone and using the 
Ge form factor given in [2^, we find a minimum upper 
limit of 2.7 X 10-38 cm^ (1.8 x lO'^^cm^) at 90% C.L. 
for this data set (combined Soudan data). 

CDMS has maintained high dark matter discovery po- 
tential by limiting expected backgrounds to less than 
one event in the signal region. These results from our 
Soudan measurements set the best WIMP sensitivity for 
spin-independent WIMP-nucleon interactions over a wide 
range of WIMP masses. Our new limits cut significantly 
into previously unexplored regions of the central param- 
eter space predicted by supersymmetry. 
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FIG. 4: Spin-independent WIMP-nucleon cross-section up- 
per limits (90% C.L.) versus WIMP mass. The upper curve 
(dash-dot) is the result of a re-analysis [l7j of our previously 
published data. The upper solid line is the limit from this 
work. The combined CDMS limit (lower solid line) has the 
same minimum cross-section as XENONIO 19] (dashed) re- 
ports, but has more sensitivity at higher masses. Parame- 
ter ranges expected from supersymmetric models described 
in [531 (grey) and [2l[ are shown (95% and 68% confidence 
levels in green and blue, respectively). Data courtesy of [2^ . 
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